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Abstract 

A set of quantum error correcting codes based on classical Reed-Muller 
codes is described. The codes have parameters [[n, k, d]] = [[2 r , 2 r — C(r, t) — 
2Ei=oC(M), 2* + 2* 



it-n 



The study of quantum information is currently stimulating much interest. Most 
of the basic concepts of classical information theory have counterparts in quantum 
information theory, and among these is the idea of an error correcting code. An error 
correcting code is a means of storing information (whether quantum or classical) in a 
set of bits (ie either qubits or classical bits) in such a way that the information can be 
extracted even after a subset of the bits has changed in an unknown way. Such codes 
are a fundamental part of the study of classical information channels. The possibility 
of quantum error correction was only recently discovered |lj, ||| . Importantly, it was 
shown that efficient quantum codes exist for arbitrarily large amounts of quantum 
information || |]]. The word 'efficient' refers to the fact that the rate k/n of the 
code need not fall off as n increases, for a given ratio d/n, where d is the minimum 
distance of the code. This and other features makes quantum error correction the 
best prospect for enabling quantum information to be transmitted or stored with a 
small amount of error, and consequently the best prospect for controlling noise in a 
quantum information processor. 

The subject of quantum error correction may be considered to have two distinct 
parts. The first part is to show how to apply error correction in a physical situa- 
tion, and the second is to find good quantum error correcting codes. This paper is 
concerned with the second part, that of finding codes. 



1 



Following 0, I will use the notation [[n, k,d}} to refer to a quantum error cor- 
recting code for n qubits having 2 k codewords and minimum distance d (previously 
I used {n, k, d} |l I]). Such a code enables the quantum information to be restored 
after any set of up to [(d — 1)/2J qubits has undergone errors. In addition, when d 
is even, d/2 errors can be detected. Attention will be limited to the 'worst case' that 
any defecting qubit (ie any qubit undergoing an unknown interaction) might change 
state in a completely unknown way, so all the error processes X, Z and Y = XZ 
must be correctable 0, ||, [7], [12j , where 



X 




(1) 



The first known general code construction used a pair of classical codes C\ = 
[n,ki,d], C 2 = [n,k2,d] with C 2 < C\ to produce a quantum [[n,ki + k 2 — n,d}} 
codeQ 0, |]. Further individual quantum codes were subsequently found || |, |H], PH 
which are more efficient than those given by this construction. Gottesman found 
an infinite set of optimal single-error correcting quantum codes, with parameters 
[[2 r ,2 r — r — 2,3]]. In this paper I will present a set of quantum codes of which 
Gottesman's is a subset. They are obtained by combining classical Reed-Muller 
codes 



13 , and have parameters 



[[n, k, 



t-i 



2 r , 2 r -C(r,t) -2^C(r,i), 2* + 2*" 1 

i=0 



(2) 



where C(r,t) = r\/t\(r — t)\. These parameters are tabulated in table 1 for small 
values of r and t. 



To show how the codes are derived, we will use the results and notation of 



Calderbank et al |10| who showed how to reduce the quantum coding problem to 
one of orthogonal geometry. A quantum code for n qubits is specified by its generator 
matrix which has the general form 



G X \G Z 



(3) 



where G x and G z generate n-bit binary vector spaces. The rows of G x and G z have 
length n, and the number of rows is n + k. The minimum distance of the code is the 
minimum weight of a non-zero generated code word, where the weight is the number 
of non-zero bit locations. A bit location is non-zero if it is non-zero in either the left 
hand (X) part or the right hand (Z) part. In other words, if a code word is written 

1 The 'less than' sign is used here to mean that one linear vector space is a subspace of the 
other. 
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where u x and u z are n-bit strings, then the weight is the Hamming weight 



\u x \u z 

of the bitwise OR of u x with u z 



To qualify as an error correcting code, the quantum code must satisfy a property 
which is best specified in terms its stabiliser 



n 



The stabiliser was discussed in |T2| , it is related to the generator by |]TD| 

0. 



H x ■ G z 



(4) 



(5) 



where the arithmatic is over a binary field (ie ■ is bitwise and and + is bitwise xor). 
From this relation, it is clear that the stabiliser and generator are the quantities 
corresponding to the parity check and generator matrices for a classical code (relation 
(H) states that 7i may be obtained from Q by swapping the X and Z parts, and 
extracting the dual of the resulting (n+k) x In binary matrix). The further property 
which a quantum code must satisfy is IfLOfl 



0. 



(6) 



The encoding method of || |J using pairs of classical codes leads to a generator 
and stabilisor of the form 




n 





(7) 



where G\ and G2 generate the classical codes C\ and C%. It is seen that the relation 
(ID is the generalisation of the 'dual' condition < C\ which such codes must 
satisfy. Previously I suggested extending these codes by adding further rows D 
to Gi such that G\ and D together generate a classical code of smaller minimum 
distance than G\ alone, and then introducing sign changes among the code vectors 



in a systematic way |11]. This method corresponds to a generator of the form 



G l 








G 2 


D x 


D z 



(8) 



11] becomes D x here. This method is adopted in what 



where the matrix D of 
follows. 



The following procedure will generate a quantum code having parameters given 
in equation (^|). Explicit examples are given after the general construction. In 
equation take 



G 1 = G 2 = 2 r , k RM (t,r), 2 t+i Reed-Muller code, 



(9) 
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where 

k RM {t,r) = 2 r -J2C(r,i). (10) 

Take D x such that G\ and D x together generate the Reed-Muller code of distance 
2*. This classical code has size /crm(^ — l, r ) so the total number of rows in the 
quantum generator (equation (||)) is 

n + k = k RM (t,r) + k RM (t- l,r) (11) 

which gives the size k of the quantum code in equation (^). If D z were zero, then 
the generated quantum code would have minimum distance 2 t . We now construct 
D z so as to increase this minimum distance to 2* + 2* _1 . Let D x refer to the z'th row 
of D x , and let the rows be numbered 1 to m where m = k RM (t — 1, r) — k RM (t, r). 
Form D z such that 

D\ = D x +1 for 1 < i < m , . 

where the action of L t is to rotate the bit string left by t places (eg L 2 (00010011) = 
01001100). This completes the construction of the generator Q, which completely 
specifies the quantum code. 

The proof that these codes have minimum distance 2* + 2 l ~ l for all r is by 
induction on t, using the fact that this was proved for t = 1 previously. The 
induction on t (ie the fact that 2* + 2 t ~ 1 is the correct formula for d) is obvious from 
the construction of the codes. An alternative proof should be possible by induction 
on r, following a similar method to a proof for classical Reed-Muller codes |T3[ 



though I have not yet established this explicitly. To prove that the codes satisfy the 
'self dual' condition (|6p, notice that the stabiliser H of a code [[n, k, d]) constructed 
as above is equal to the generator of a [[n, —k, d']] code constructed by the same 
method (where d' = 2d, (4d) for odd, (even) r respectively, see table 1). In other 
words, the stabiliser matrix is also built out of classical Reed-Muller matrices by 
the same recipe. Because the classical Reed-Muller codes contain their duals, it is 
easy to show that whenever the number of rows in the stabiliser is less than n (ie 
k > 0), H x ■ if J = and therefore equation ([]) is satisfied (illustrative examples will 
be given shortly). When the stabiliser or generator matrix has n rows, on the other 
hand, one finds D x ■ D T Z ^ which, combined with the fact that the rest of H x ■ if J 
is zero, means that H x ■ ifj cannot be symetric and therefore @ is not satisfied. In 
summary, the method works for k > 0, and the generators for k < are stabilisers 
of k > codes. 

Some examples will clarify the construction. 
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First take the simplest case r = 2, t = 1. This produces the following generator 



/ 1111 
0000 
0011 

V oioi 



0000 \ 

1111 

0101 
0110 / 



(13) 



which does not satisfy (|6]), but is significant in that it 'seeds' the rest of the con- 
struction method. 

Next, with r = 3, t = 1 we obtain the [[8, 3, 3]] code first derived independently 
by several authors [|1^, [H], [T| : 



£[[8,3,3]] 



/ 11111111 
00001111 
00110011 
01010101 
00000000 
00000000 
00000000 
00000000 
00000011 
00000101 
V 00010001 



00000000 \ 
00000000 
00000000 
00000000 

11111111 

00001111 
00110011 
01010101 
00000101 
00010001 
00000110 / 



(14) 



The whole set of codes with t = 1 are those derived previously by Gottesman [12 . 
The present discussion constitutes a particularly straightforward way of constructing 
them. The stabiliser of such a code is 



n 



[[8,3,3]] 



/ 11111111 

00000000 
00001111 
00110011 

V oioioioi 



00000000 \ 

11111111 

00110011 
01010101 
00111100 J 



(15) 



To see more clearly why the method fails for k = 0, consider the following 
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stabiliser for t = 2, r = 4: 



:2,r=4 



/ 1111111111111111 

0000000011111111 
0000111100001111 
0011001100110011 
0101010101010101 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000001111 
0000000000110011 
0000000001010101 
0000001100000011 
0000010100000101 
V 0001000100010001 



0000000000000000 \ 
0000000000000000 
0000000000000000 
0000000000000000 
0000000000000000 

1111111111111111 

0000000011111111 
0000111100001111 
0011001100110011 
0101010101010101 
0000000000110011 
0000000001010101 
0000001100000011 
0000010100000101 
0001000100010001 
0000000000111100 / 



(16) 



Comparing with equation (|I5|) , one sees that the problem is that the ones are too 
spaced out in the final 3 rows of H x for this part to satisfy the parity checks of the 
previous 3 rows. This happens whenever the relevent classical Reed-Muller code is 
too small to contain its dual, which is the case for k < in our construction. 



Next, 
table in 



the stabiliser for t = 2, r = 5 will be shown. It is close to optimal (see 
D), and is the simplest previously unknown code obtained in this paper: 



Ti 



[[32,10,6]] 



(17) 



£[32,16,8] 


0000000000000000000000001 1 1 1 1 1 1 1 
00000000000000000000111100001111 
00000000000000000011001100110011 
00000000000000000101010101010101 
00000000000011110000000000001111 
00000000001100110000000000110011 
00000000010101010000000001010101 
00000011000000110000001100000011 
00000101000001010000010100000101 
V 00010001000100010001000100010001 



\ 

G[32, 16,8] 
00000000000000000000111100001111 
00000000000000000011001100110011 
00000000000000000101010101010101 
00000000000011110000000000001111 
00000000001100110000000000110011 
00000000010101010000000001010101 
00000011000000110000001100000011 
00000101000001010000010100000101 
00010001000100010001000100010001 
00000000000000000000111111110000 J 

(18) 
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Before concluding, let us note that many other good quantum codes can be 
obtained either directly or indirectly by the construction shown in equation (H). I 
have thus been able to obtain many of the codes tabulated by Calderbank et. al. 
[§]. An example of a code obtained indirectly is the [[5,1,3]] perfect code, which 
is obtained by deleting any bit from the following [[6,0,4]] code constructed using 
equation (||): 

/ 001111 000000 \ 
110011 000000 
000000 001111 
000000 110011 

linn oioioi 

V 010101 100101 j 



^[[6,0,4]] 



(19) 



I will conclude with a few comments on the usefulness of the quantum Reed- 
Muller codes. Like their classical counterparts, these codes are easily encoded and 
decoded. This may be especially important for applications in a quantum computer 
where efficient correction is neccessary. The codes are optimal (have maximal k for 
given n and d) for t = 1 [12|, and close to optimal for small t and r (comparing with 
the table of codes in [El, it is seen that [[32, 10, 6]] is close to optimal). They are far 
from optimal for large t and r, but should form a good starting point for deriving 
other codes. The code size falls to zero at t — r/2 with even r, giving parameters 
[[2 r , 0, 2 r l 2 + 2 r ' //2 ~ 1 ]], ie d ~ n 1 ^ 2 for n ^> 1, which means that in the limit of large 
r the small size codes are no better than those obtained by the simpler method 
of equation (|7|) applied to the self-dual Reed-Muller codes, leading to parameters 
[[2 r - 1 ,0,2 p / 3 ]]. 
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94 


-±o 


A 
~± 

8 


n 

3 


-3 








16 


10 





-10 






n 32 


25 


10 


-10 


-25 




64 


56 


35 





-35 


-56 


128 


119 


91 


35 


-35 


-91 


256 


246 


210 


126 





-126 


512 


501 


456 


336 


126 


-126 


1024 


1012 


957 


792 


462 






Table 1. Parameters [[n, k, d\) of the quantum Reed-Muller codes, given by equa- 
tion ©, for 2 < r < 10, 1 < t < 5. 
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